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We investigate whether a mass scale for elementary particles can be derived from interactions 
of particles with the distant matter in the Universe, the mechanism of the interaction being the 
classical vector potential, propagating in a space of negative curvature. A possible context for such 
a mass scale is conformal gravity. This theory may prove to be renormalizable, since all coupling 
constants are dimensionless; conversely, however, there is no coupling constant analogous to the 
conventional G to provide a starting point for a mass scale calculation. We obtain the equations 
for propagation of the vector potential of a charged particle moving in a plasma in a curved space. 
We then show that distant matter will contribute to (A^A^), and that this non-thermal part will 
eventually dominate the ordinary thermal part. At this point a symmetry breaking transition of the 
Coleman- Weinberg type is possible, and particle masses can be generated with m 2 « {A^A^). 
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I. INTRODUCTION 

In this paper we will try to implement a possibility 
suggested in a well-known text on quantum field the- 
ory [1]. After surveying the difficulties faced by current 
theories of the mass scale of elementary particles, the 
authors write: ". . . it may be that the overall scale of 
energy-momentum is genuinely ambiguous and is set by 
a cosmological boundary condition." The difficulties the 
authors are concerned about begin to appear in the usual 
description of the origin of the mass scale: 

1. The Einstein theory of gravity, based on an action 
that is linear in the curvature tensor, is the correct 
one. 

2. The coupling constants of the fundamental interac- 
tions are all dimensionless except for the gravita- 
tional constant, G. As a result, all interactions are 
renormalizable except for gravity. 

3. G can be combined with other constants to define 
a mass, the Planck mass, m p i anc k- 

4. Because gravity is not renormalizable, the masses 
of particles will be dominated by high-energy gravi- 
tational self-interactions, and the masses should be 
of the same order as the Planck mass. 

This brings us at once to the "hierarchy problem" : why 
is the ratio m pr oton/n p i arlc k so small, about 10 -19 ? 

In this paper we present a very different model of the 
origin of the mass scale, with these main features: 

1. A gravitational theory that is renormalizable, and 
has dimensionless coupling constants, will ulti- 
mately come to be accepted in place of the Einstein 
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theory. A possible example of such a theory is con- 
formal gravity, which is based on an action that is 
quadratic in the curvature tensor (see [2] , and other 
papers cited there). It is too early to say whether 
conformal gravity will be able to describe recent 
cosmological observations, but we will assume that 
whatever theory is ultimately accepted will have 
the basic characteristics introduced here. In most 
of this paper the precise form of the gravitational 
field equations will be irrelevant. We will only need 
these equations in appendix B, where we consider 
Mannheim's model [2] of conformal gravity. 

2. In a renormalizable theory the mass scale may not 
be generated by high-energy virtual processes at 
all, and we suggest that it actually arises from in- 
teraction with distant matter in the Universe. 

3. The agent of this interaction is the familiar clas- 
sical vector potential of the electromagnetic field, 
propagating in a space of negative curvature. 

4. The influence of distant matter rises steadily from 
the birth of the Universe until a symmetry-breaking 
transition takes place, possibly of the Coleman- 
Weinberg type [3], when masses as we know them 
appear. We use the letters CW, in text and sub- 
scripts, to refer to this transition. 

We start from the simple, and at first sight pointless, 
observation that the kinetic-energy term of the Klein- 
Gordon equation for a charged scalar field in Minkowski 
space, 4>*{\> — qA) 2 cj), when expanded, gives the term 
(f)*q 2 A 2 (f>. This has the same sign as the mass term 
4>*m 2 4>, suggesting that in some circumstances, such as 
the radiation era of the early Universe, (q 2 A 2 ) might play 
the role of m 2 . In the familiar world of of Minkowski 
space this will not happen, because (A 2 ) depends sim- 
ply on the local temperature. When pulses propagate in 
a curved space, however, the potentials (though not the 
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fields) leave a tail, as will be shown in detail in section IX. 
As a result, distant matter will generate a non-thermal 
component in (A 2 ). 

Our model will use the preferred geometry for con- 
formal cosmology [4], a Friedmann- Robertson- Walker 
(FRW) Universe with negative curvature. The current 
consensus among cosmologists is that space is flat, at 
least at the present time. But we must remember that 
the common inference that space is flat rests on the usual 
assumptions about gravity and the generation of mass. 
In this paper these assumptions are in abeyance, so the 
curvature of space remains an open question. 

The early Universe is filled with radiation-dominated 
plasma, in which each charged particle is accompanied by 
a screening cloud. It is by no means obvious how elec- 
tromagnetic interactions could extend over cosmological 
distances, rather than declining exponentially over dis- 
tances of the order of a Debye length. In section IV we 
explain how this happens, and describe our model of the 
sources of these long-range fields and potentials. 

In subsequent sections we study the propagation of 
electromagnetic waves in curved space. Here we will not 
need any gravitational field equations, but will simply 
use the well-known equations for propagation of a classi- 
cal electromagnetic field [5]. However, we are assuming 
that the vector potential can in some circumstances act 
as a direct agent, rather than through its derivatives, so 
the gauge cannot be freely chosen and we will have to 
pay careful attention to the gauge condition. We discuss 
gauge invariance and the gauge condition later, in section 
VIII. 

In sections IX and X we give the simplest derivation of 
the mass scale, assuming that the long-range fields from 
the plasma are not thermalized, even when propagating 
in a curved space. 

Sections XI through XV can be omitted on first read- 
ing. Here we assume the long-range fields are slowly ther- 
malized, and verify that even in this situation the mass 
scale will be established. 

Section XVI gives our conclusions. 



II. NOTATION 



The expansion parameter, R(t), has the dimension of 
length, and can be thought of as the radius of the Uni- 
verse at time t. 

The conformal time, 77, is related to the ordinary time t 
by di] = dt/R{t). In terms of ij, x, 9, 4> (all dimensionless) 
we have: 

ds 2 = R 2 (r))(- dry 2 + d X 2 + sinh 2 x d<9 2 

+ sinh 2 xsin 2 (9d0 2 ) , (4) 
9^u = i? 2 (r7)diag(-l,l,sinh 2 x,sinh 2 xsin 2 6>) , (5) 
V5 = R 4 (v) sinh 2 x sin 9. (6) 
Maxwell's equations in free space are: 

(VgF af3 ),p = 0, (7) 

-Fa/3, 7 + Fp*f,a + F~/u,f3 = 0. (8) 

Potentials and the gauge condition will be introduced 
later. 



III. THE MODEL 

We will assume the Universe has the form of a bub- 
ble that condenses out of a metastable exterior vacuum 
[7, 8] . The surface of the bubble expands at the speed of 
light, and the interior of the bubble forms a FRW space 
with negative curvature. We are not concerned with in- 
flation, but will simply assume that the material content 
of the Universe appears as radiation at the surface of the 
bubble. 

We restrict the material content of our model to the 
following: 

1. A charged, massless, scalar field, <fr, with conformal 
weight —1. 

2. The usual electromagnetic fields and potentials. 

3. In appendix B we implicitly assume, following 
Mannheim [2], an additional neutral scalar field, S, 
which from the beginning has a very large, constant 
value, So, but will play no further role. 



Units are chosen so that c = ft = 1. Our sign conven- 
tions are those of Weinberg [6] , so the metric signature 
is (- + ++) and g = -Dct (,g M „). 

For the early Universe, assumed spatially homoge- 
neous, we use a FRW metric, with coordinates t, x, 9, (f>, 
labeled 0, 1, 2, 3: 

ds 2 = -dt 2 + R 2 (t) (dx 2 + sinh 2 xd6» 2 

+ sinh 2 xsin 2 6»d0 2 ) , (1) 

9llv = diag(-l,i? 2 (t),i? 2 (t)sinh 2 x, 

i? 2 (t)sinh 2 xsin 2 <9) , (2) 
y/g = i? 3 (t)sinh 2 x sin 6>. (3) 



IV. THE SCREENING CLOUD 

We will be concerned with the propagation of the vec- 
tor potential generated by moving charges in the plasma 
of the early Universe. The charges will have a thermal 
mass due to their interaction with the photons, and the 
plasma will be relativistic. For simplicity, however, we 
will carry over some results from non-relativistic plas- 
mas. 

A stationary charged particle in the plasma will be 
surrounded by a screening cloud that reduces the field 
exponentially on a scale known as the Debye length ([9], 
chapter 2). Surprisingly however, starting in the 1950's 
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several authors [10-12] discovered that a test charge mov- 
ing uniformly is not exponentially screened, but gener- 
ates the field of a quadrupole in a collisionless plasma. 
The restriction to a collisionless plasma was lifted in later 
papers [13, 14], where it was shown that the far field of 
a moving test charge is that of a dipole with strength of 
order qr t V, where q is the charge, r t the collision time as 
measured by t, and V the velocity. 

A charge that is part of the plasma, rather than being 
a test charge constrained to move in a certain way, will 
itself experience collisions. We can picture such a particle 
as describing Brownian motion, with its screening cloud 
sometimes closer, sometimes farther, but never fully es- 
tablished. We will model the most important aspect of 
a single step of this process by setting up a local set of 
Cartesian coordinates and placing a stationary screening 
charge — q at the origin. A charge q is initially at rest 
on the z axis at z = —Vr t , then begins to move with 
uniform velocity V along the z axis, starting at t = —r t 
and ending at t = Tt, when it again comes to rest. Note 
that this choice of origin for t differs from that of the cos- 
mic time used in section III; n, defined below, will differ 
in a similar way. We reconcile these differences later, in 
section X. 

We are now going to take the result for a uniformly 
moving test charge and use it to calculate the field around 
our model charge that is subject to Brownian motion. We 
justify this step as follows. Imagine that the test charge, 
instead of moving uniformly from an indefinite time in 
the past, is actually stationary until time t = 0, and then 
starts moving uniformly. There will be transient fields, 
but after a few collision times the final field of a moving 
dipole will be established. During the transient period, 
the test charge will be moving away from its screening 
charge, which only gradually picks up speed and trails 
along behind. This is similar to the motion of the charge 
in our model. Since the transient fields form a bridge be- 
tween the initial and final fields they must have a similar 
form to the field of the test charge, and extend compara- 
bly far. 

During the period that the charge is moving (the 
pulse), the dipole moment will be D(t) = qVt. Before 
and after the pulse the dipole moment remains at a con- 
stant value, but this is of little interest because a constant 
dipole moment generates no magnetic field or vector po- 
tential. It is well known that Maxwell's equations sepa- 
rate in conformal time, so we will write the dipole mo- 
ment in terms of n. Define t = Tt/R s (t) (dimensionlcss) , 
so that the pulse extends from rj = — r to n = r. The 
subscript 's' on R s indicates the source time. R s (t) can 
be treated as constant during the pulse, and the dipole 
moment can be expressed as a Fourier integral: 

D(r]) = 77- / D(n) exp(— inrj) dn, (9) 
2tt J 

2iqVR s sin(n-r) 



D(n) 



(10) 



of the local system used in this section will be parallel to 
the polar axis of the polar coordinates used in the bulk 
of the paper. 



V. DIPOLE FIELDS 

We will need only the three fields F w , F 2 a and F i2 . 
Since Maxwell's equations separate in conformal time the 
elementary solutions can be written 

F w = .fio(x,^,0)exp(-mr?), (11) 
F20 = /2o(x,M)exp(-mr7), (12) 
F12 = /i 2 (x,M)exp(-mr7). (13) 

For dipole fields we try the following forms: 



ho = .h(x)Pi(0), 

/20 = f2(x)Ni(6), 
fi2 = Mx)Ni(6), 



(14) 
(15) 
(16) 



with angular functions P\{9) = cos 9 and Ni(6) = 
dPi(0)/d0 = -sin9. 

Maxwell's equations now give (with a prime meaning 
d/d X ): 



msinh x/i - 2/ 3 = 0, 
W/2-/3 = 0, 
fi-fi+infa = 0. 

A. Dipole magnetic field 



(17) 
(18) 
(19) 



From (17), (18) and (19) we obtain the equation for / 3 
alone 



d 2 / 3 
dx 2 



sinh 2 x 



h = 0. 



(20) 



This is the analog of equation (16) of Mashhoon [15], 
which he derived for a space of positive curvature. We 
solve (20) by first expressing the solution in terms of 
hypergeometric functions [16]. We then show that these 
particular hypergeometric functions can themselves be 
expressed in closed form in terms of simpler functions. 
The indicial equation is p(p — 1) — 2 = 0, so p = 2 or 
p = — 1. Write the regular and singular solutions as 



/ 3 , reg = sinh 2 xP 3 ,rog, 
/3,sin g = sinh _1 xP 3 , s ing- 



(21) 
(22) 



We estimate the value of r t in section X. The z axis 



It is then straightforward to show that, in terms of the 
variable /j = (1 — coshx)/2, 

F 3;reg = F(2 + in,2-in;5/2;n), (23) 
Ps.sing = F(-l + in, -I -in; -1/2; n). (24) 

Explicit closed forms for these hypergeometric func- 
tions are given in appendix A. / 3 is constructed from 
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that combination of fa^eg and facing that is proportional 
to exp(inx), because when combined with exp(— inr/) this 
represents outgoing waves: 

Ci(n)exp(mx) , . . . , . 
sinhx (coshx-msinhx), (25) 

where Ci(n) is a normalizing factor to be determined. 

It is convenient at this point to express fa in terms of 
u = tanh(x/2): 

CAri) exp(inx) / „ 9 \ ,„„s 

fa = 1 ' 2u ' (l-2inu + u 2 ). (26) 



B. Dipole electric fields 

From (17) and (25) we derive the equation for the ra- 
dial electric dipole field: 

/i = - tCl{n) eX ^T )(1 = (1 - 2inu + u 2 ) . (27) 

Similarly, from (18) and (26) we derive the equation 
for the transverse electric field: 



h = 



iC\{n) exp(mx) 
4nu 2 

x [(1 - u 2 ) 2 - 2inu{\ + u 2 ) - An 2 u 2 ] . (28) 



C. Normalizing the fields 

From (9) and (10) we can get the near-field expression 
for the radial component of E, and by comparison with 
(27) arrive at the form of the normalizing factor C\ (n) . 
On the polar axis, at small distances, the radial E field 
at frequency n is 



Er 



d f D{n) exp(— inn) 
dr \ r 2 
2D(n) exp(~inr]) 



(29) 



where r = R s x- 

In the same limit (small \), (27) gives /i(x) = 
— 2iCi(n)/(nx 3 ), and so, on the axis, 

Fiotn, X) - -2zC 1 (n)i? s 3 cxpi-m^/inr 3 ). (30) 

We now transform to a local Minkowski frame with 
coordinates t, r: 

= -2iG\(n)R s aq>{-inr))/{nr 3 ). (31) 

Fio is just the conventional radial electric field as given 
in (29), so using (10): 



Ci(n) = 



— 2gVsin(nr) 



(32) 



VI. PROPAGATION OF THE MAGNETIC 
FIELD 

(26), (32), (16) and (13) give the Fourier transform of 
the magnetic field: 



F12 — 



—qVsin(nT)Ni(9) exp[in(x — f])] 



nu 

2\ 



x (1 - 2inu + u 2 ) . 



(33) 



We transform back into x, V space by dividing by 2ir 
and integrating over n along the real axis. For \ > rj + r 
both exponentials in sin(nr) allow us to close in the UHP, 
using the contour of figure 1. sin(nr)/n is regular at 
n = 0, so we get zero, as required by causality. Similarly, 
when x < t]—t we can close in the LHP, using a contour 
that is the inverse of figure 1, and we again get zero. 




FIG. 1: Typical contour for Fourier synthesis of magnetic and 
electric fields. 

For Tj — t < x < ?? + Twe must write sin(nr) = 
(exp(mr) — exp(— inr))/2i, divide the integrand into two 
pieces, and close the first integral in the UHP and the 
second in the LHP. Combining these two integrals we get 



F12 — 



qVsm6(l + u 2 ) 
2~u ' 



(34) 



The propagation of F\i is shown in figure 2. An un- 
expected feature of the pulse is that the amplitude does 
not tend to zero for large x, but to a constant asymptotic 
value. The conventional H field, docs, of course, tend to 
zero, in fact exponentially, because in a local Lorentz 
frame we have F\i = F\i /{R 2 sinh x) , the subscript 'p' 
on i? p indicating the point of observation. 



VII. PROPAGATION OF THE ELECTRIC 
FIELD 

The propagation of the electric field can be displayed 
in a similar way, with the contour determined by the 
boundary conditions for r\ < — t, before the pulse begins. 
Referring to figure 1, we have to restrict our integral to 
the portions A-B and C-D, excluding the semicircle B- 
C. The integral is interpreted in a principal value sense. 
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To fix the potentials uniquely we need a gauge con- 
dition. The one often suggested is the Lorenz condition 
[5, 17]: 



(42) 



But this is not acceptable here, because we want a 
theory that is conformally invariant, and the Lorenz con- 
dition is not (the conformal weights of y/g, and 
are 4, —2 and 0, respectively). We propose in this paper 
to use the modified condition 



o. 



(43) 



This is conformally invariant since <f> has conformal 
weight —1. In applying (43) we set <j)*<j) equal to its ex- 
pectation value, which is proportional to T 2 , and hence 
to i?~ 2 . All factors of R now disappear from (43), and 
we obtain 



FIG. 2: Propagation of the magnetic field, F12: snapshots of 
the pulse for various values of conformal time, r\. For each r\, 
the pulse is confined to the range r) — t < \ < V + T '> within 
that range we plot (1 + u 2 )/u. To render the pulse visible on 
the graph, we have arbitrarily set r = 0.5, far larger than it 
is in practice. 



We now find there are electric fields both before and after 
the main pulse. This is to be expected, because in our 
model the source has a static dipole moment before the 
pulse begins, and an opposite moment after it has ended. 



in sinh 2 xho + -7- [sinh 2 \hi\ - 2h 2 = 0. (44) 

A 

The potentials are still not defined uniquely. We are 
permitted to make a restricted gauge transformation, 
with a gauge function ip that satisfies 



0. 



(45) 



We will deal with this remaining ambiguity below, in 
subsection VIII B. 



VIII. DIPOLE POTENTIALS 



A. Scalar potential 



The potentials are defined by the following equations: 



A) = h o (x)Pi(0)exp(-inri), 
A 1 = fci(x)Pi(0)exp(-inf7), 
A 2 = h 2 (x)N 1 (6)exp{~inr)), 

dA v dA^ 

From (38) we derive: 

h>2-h! = h, 
Hq + inh\ = /1, 
h Q + inh 2 = f 2 . 



(35) 
(36) 
(37) 

(38) 



(39) 
(40) 
(41) 



(39), (40) and (41) are not independent. If we differen- 
tiate (41), subtract it from (40) and use (19), we obtain 
(39). These three equations therefore do not suffice to 
define the potentials uniquely; this is to be expected, be- 
cause A^ is only defined up to a scalar gauge function, 
V>, so that A^ = A^ + dip/dx^ defines the same fields as 
A,,. 



By combining (17), (18), (40), (41) and (44) we can 
derive an equation for ho alone: 

d 2 h „ , dh ( 2 2 \, n 

" +2 coth X -7 1 + [n 2 -— — Wi = 0. 46 
dx dx V sinh xJ 

This is the equation for hyperbolic spherical functions 
[8, 18]. We will solve it by means of hypergeomctric 
functions, as for the fields. The indicial equation gives 
p = 1 or p = —2, and the two solutions can be shown to 
be: 

fto.reg = smhxF(2 + ia,2-ia;5/2;n), (47) 
Vsing = sinh~ 2 X F(-l + ia, -I -ia; -1/2; (i). (48) 



where [i = (1 — coshx)/2 and a = \Jn 2 — 1. The hyper- 
geometric functions are the same ones we encountered 
with /3, but the powers of sinhx appearing in front of 
them are different. 

A will be constructed from that combination of the 
two solutions that represents outgoing waves, at least for 
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large n, where a is real. With u = tanh(x/2), as before: 

A a = Pi(0)/io(x)exp(-mr?) 

= C 2 {n)Pi (9) exp[i(«x - nrj)} 

(1 - u 2 ) (l - 2iau + u 2 ) , 
x [ - l±- a >-, (49) 

where C2(n) is a normalization factor that depends on 
the driving function. C 2 (n) is, of course, proportional to 
the normalization factor C\{n) defined in (32). For small 
X, (49), (41) and (28) give: 

*Ci(n) 



C 2 (n) 



An 



(50) 



B. Restricted gauge transformations 



For restricted gauge transformations, a suitable gauge 



function of dipole form is ip 
where h satisfies 

d 2 h n , dh 
—- + 2cothx-^- - 

ox ox 



Kx)Pi{0) exp(-mr?), 



sinh z 



X 



h = 0. 



(51) 



This is the same equation that is satisfied by the scalar 
potential. If we include tp in our definitions of the poten- 
tials, ho, for example, becomes ho = ho — inh. But ho 
permits no such addition; it is already completely speci- 
fied by causality and normalization. 



C. Vector potential 

The interesting component of the vector potential is 
the transverse one, Ag, or in our notation A 2 . (41) gives 



hn = 



-if 2 iho 
n n 



(52) 



This can be used to get hi, and from that A 2l since f 2 
and ho are known from (28) and (49). 

Equation (52) shows that A 2 is the sum of two parts, 
one involving exp(inx), the other involving exp(iax): 

M = w {Af+A) ' (53) 

N = e m *[{l-u 2 ) 2 -2inu{l + u 2 )-An 2 u 2 } , (54) 



A = e lax (l-u 2 ) [-l + 2iau- 



(55) 



D. Gauge invariance in the end 



We have set up the scalar and vector potentials accord- 
ing to a "classical prescription" : treat the charges in the 
plasma as classical particles, and use the simple criteria 
of causality and conformal invariance to fix the poten- 
tials uniquely. Once we have done this, we can represent 
the charged particles as quantum fields. General gauge 
transformations are then permitted, provided the phases 
of the wave functions of charged particles are simultane- 
ously transformed in the usual way. 



IX. PROPAGATION OF THE VECTOR 
POTENTIAL 



The Fourier synthesis of A 2 is carried out, as for the 
magnetic field, by dividing by 2n and then integrating 
over n, along a contour chosen to respect causality. There 
is no vector potential before the pulse begins, so the cor- 
rect contour is a line parallel to the real axis and slightly 
above it. The integral must be taken along the whole 
path A-D in figure 1, including the semicircle B-C. This 
ensures that for x > V + T the contour can be closed in 
the UHP and the integral will be zero. There arc two 
regions of interest, i]—T<x<il + T (the "main pulse"), 
and x < V ~ T (the "tail" ) . 

We get a non-zero result only for that part of the inte- 
gral that involves a contour that is closed in the LHP. For 
the part of equation (53) that is proportional to exp(inx) 
we can shrink this contour to a small circle about the 
point n = 0, and we just have to find the residue there. 
For the part that is proportional to exp(iax) we have to 
remember the branch points at n = ±1, so our contour 
can only be shrunk to the form shown in figure 3. 



FIG. 3: Integration contour for A2. The contour encloses the 
points n = — 1 and n = +1; it is described in a clockwise 
sense. 



A 2 is the sum of contributions from small circles (or 
semicircles) around n = (the pole terms) and the in- 
tegrals along the cut from n = — 1 to n = 1. The pole 
terms have simple analytic expressions, but the integrals 
must be evaluated numerically for each (77, x) pair. The 
integrations are straightforward, and the propagation of 
A 2 is shown in figure 4. 

An important difference between figures 2 and 4 is that 
in the latter the pulse has a non-zero tail for x < V — t. 
This feature of the propagation of potentials in a curved 
space has been noted before [19, 20]. For dipole prop- 
agation, as here, the tail rises linearly from x = 0, and 
approaches a constant value for large x- This asymptotic 
value is proportional to r. In ordinary Minkowski space, 
which corresponds to the limit x - * 0, A has no tail. 

We note that the following simple function, with a = 
0.735, gives an adequate fit for A\ for large rj: 



(2qVr sin6») tanh 2 (ax). 



(56) 



7 



FIG. 4: Propagation of the transverse component of the vec- 
tor potential, A2 : snapshots of the pulse and the tail for a se- 
ries of values of 77. As with the magnetic field, we set r = 0.5 
to make the pulse visible on the graph. 



MASS GENERATION: 
CALCULATION 



FIRST 



Let us temporarily set aside all questions of gauge in- 
variance, and simply regard (A^A^) as being composed 
of a thermal part, (A^A^th, which will be of order T 2 , 
and a non-thermal part, (A^A^nt. Moving particles in 
the distant plasma will generate pulses of A^, each of 
which consists of a "main pulse" and a "tail" . The main 
pulse contains electric and magnetic fields, and will con- 
tribute to {A^A^th- The tail, however, is a pure gauge 
potential that generates no fields; it will contribute to 
<AM„) nt . 

The development of (A^A^)^ can be visualized as fol- 
lows. As each pulse passes the observation point it leaves 
a memory in the form of the tail. These tails will not 
cancel but will add according to the theory of random 
flights [21]. The scalar potential, of course, will remain 
close to zero; it is only the vector potential that accu- 
mulates these additions. Consequently, from now on we 
will write (A^A^th — (A 2 ) t h, and similarly for the non- 
thermal part. 

We are concerned with the ratio 



p=(A 2 ) nt /(A 2 ) th . 



(57) 



Both the numerator and the denominator will decrease 
as the temperature falls, but the numerator falls more 
slowly, so the ratio will build up from zero until it be- 
comes of order unity. At this point a Coleman- Weinberg 
transition [3] can take place and normal masses will ap- 
pear. These will have magnitude m 2 w (A 2 ) nt . A 



schematic diagram illustrating this process can be found 
in Narlikar and Padmanabhan [22], figure 10.1. When 
(A 2 ) nt exceeds (A 2 ) t h a second minimum, lower than 
the one at m = 0, will develop in the effective potential. 

It is convenient at this point to change coordinates so 
the observation point is at 77 = t] p , \ = 0, and a general 
source point is at 77, v. Our formulae involve differences 
in rj and \, which are unchanged by this shift of origin. 
We will consider 77 and t as starting from zero at the time 
of minimum radius, i? m i n , and maximum temperature, 
T max . We can picture the buildup of (A 2 ) nt with the 
help of figure 5, where the coordinates are 77 (upwards) 
and X- The row of boxes at the bottom of the diagram is 
at r\ = 0; this represents the surface of the bubble. Each 
box has thickness d% and duration r. Two observation 
points are shown, Pi and P2, at times r\\ and 772- The 
lines Pi-B and P2-C represent the past light cones. All 
plasma particles within these past light cones contribute 
to (A 2 ) nt . The contributions add incoherently over times 
longer than the collision time, r, so we can imagine the 
whole diagram divided into time slices of duration t, just 
like the one shown for 77 = 0. A box at source time 773 
represents a spherical shell of volume 47ri? 3 sinh 2 \ dx- 
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X 



FIG. 5: Buildup of (A 2 ) nt . Each box represents a spherical 
shell of radius x, thickness dx, duration r. For clarity, only 
one time slice is shown; the plasma actually fills the whole 
volume from 77 = to the observation time. 

The number density of particles in this shell, for a ther- 
mal distribution, is 

u a = 0.24T 3 , (58) 
and the number of particles in the shell is then 

dN s = 47r(0.24T s 3 )^sinh 2 xdx 

= 4 7 r(0.24)4 I /4 sinh 2 xdx. (59) 

Here we have introduced, following Mannheim [2], the 
quantity Am = i? 4 T 4 , which (at our current level of 
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calculation) will be constant during the expansion, and 
therefore does not need another subscript, 's' or 'p'. 

Multiplying (59) and (56) we get the contribution to 
A\ from the particles in the shell. We convert this to 
g^A^Av by multiplying by g 22 : 



47r(0.24)< /4 sinh 2 xdx(2 g T) 2 
tanh 2 (a%) 



:(V 2 )(sin 2 6) 



R 2 sinh 2 x 



(60) 



The average of sin 2 9 over a sphere gives 2/3. For the 
average of V 2 we reason as follows: if the effective m 2 of 
the particles were purely thermal, we could write V 2 = 
1/2. But the effective mass actually increases with time 
as the ratio p from (57) increases from zero towards one. 
So a better estimate of V 2 is 1/(2 + p). 

r\ = at the surface of the bubble, and r\ — r\ v at the 
observation point. For some time slice at time r\ between 
these two limits, we get the contribution of all spherical 
shells to (A 2 ) nt by integrating from x = to the light 
cone at \ — Vp ~ V- 



M024) 3/4 



1 M 



2+p 



(2qry 



tanh 2 (ax) , , 
dx ff2 ■ (61) 



We now have to sum over all time slices of duration r. 
The sum can be converted to an integral by multiplying 
by dr]/r: 



32tt(0.24) 3/4 2 



dr/ 



1 



L 



o \2+P, 
tanh 2 (ax) 



(62) 



In section IV we defined r = T t /R s (t). For r t we will 
simply use equation (3.24) from [23]. The factor U 3 in 
the numerator will be of order (2+p)~ 3//2 ; the quantity in 
the bracket in the denominator is of order unity and will 
be omitted. We then have, in our notation, r t = to 2 (2 + 
p)~ z / 2 / (&nv s q A In A). For m 2 we will use the effective 
mass at time t. The purely thermal value is q 2 T 2 , but 
we should multiply this by (1 + p) to include the non- 
thermal part also. We have v s = 0.24T S 3 , so 



r = 



1+p 



87rln A(0.24)T s g 2 i? s (2 + p) 3 / 2 

1+p 

8tt In A(0.24) 9 2 Am /4 (2 + p)3/2 ' 



(63) 



We can now get an expression for p by using (63) in 



(62) and dividing by T 2 : 



P = 



3 In A 



% 1 



Hp 



y< 4 (2+ P )B/2 y 



L 



31nA 

4 



"■>-'» tanh 2 (a X ) 
dy 

A D2T2 

1 l p- 1 p 
1+p 

(2 + p)6/2 
1+p 



/")p 

Jo 



dry 



/ d 



dxtanh (ax) 



Vp 



drj 



31nA7 ■'(2+p) 5 /2 

x {Vp - V - tanh[a(^p - r])]/a} . (64) 

This integral equation is surely not correct in detail, 
but it does provide a general picture of the buildup of 
p from p = at r\ v = 0. We seek rj CW} the value of 
i] p for which p = 1. Use a = 0.735, as before, and for 
definiteness set In A — 10; numerical solution then gives 

Vcw « 11. 

A remarkable feature of (64) is that not only has Am 
disappeared, but also the coupling constant, q. The long- 
range character of the Coulomb interaction, however, is 
still apparent in In A. An equation of this sort might 
reasonably have been expected to yield a value of r; cw 
that was either extremely small, so the CW transition 
takes place almost immediately after the formation of 
the bubble, or extremely large, so the transition never 
takes place at all. Instead we get a value of r] cw that is 
within one or two orders of magnitude of unity. 

As an example of the application of this result, we esti- 
mate, in appendix B, parameters in the model developed 
by Mannheim [2] for conformal gravity. 

One might think that in setting up potentials by the 
"classical prescription" of subsection VIII D we would 
find that at any point a large value of (A) would de- 
velop, as well as (A 2 ) nt . But this is not so, once we treat 
the Universe as a quantum mechanical system. We must 
avoid being too specific about the state of distant mat- 
ter, which we cannot investigate directly. We should not 
assume it is in a pure state; rather, it is a mixture, and 
all states compatible with the conservation laws will be 
present in the density matrix. In many situations, as em- 
phasized by ter Haar [24] , it is a matter of taste whether 
we adopt the statistical (ensemble) or quantum mechan- 
ical interpretation of the density matrix. In the language 
of Bell [25], "or" and "and" are equally acceptable. But 
this is not the case here; the quantum mechanical inter- 
pretation is the appropriate one. In an ordinary plasma 
in Minkowski space, only the particles within the De- 
bye sphere will contribute to A, and we do not have to 
concern ourselves with mixtures. But propagation in a 
curved space obliges us to consider very large numbers of 
particles at cosmological distances, and the introduction 
of a mixture is inevitable. 

Like the distant Universe itself, (A 2 ) will be a mixture, 
and the only meaningful quantities are averages. By sym- 
metry, (A) will be zero, but (A 2 ) will not, so g 2 (A 2 ) can 
play the role of m 2 . 
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XI. INTRODUCTION OF CONDUCTIVITY 

Up to this point we have been mainly concerned with 
potentials in the tail. But we have now to look more 
closely at the electric and magnetic fields in the main 
pulse for dipole propagation. These combine to give 
a Poynting vector directed outwards. In flat space, 
this vector falls off like l/x 4 > which tends to zero even 
when integrated over the whole sphere. When we go 
to curved space, however, the Poynting vector tends 
to (qVsin9/ sinh^) 2 - The surface area of a sphere is 
4-7T-R 2 sinh 2 x, so the integral of the Poynting vector tends 
to a constant non-zero value. One effect of curvature is 
to continuously increase the energy of the plasma. 

It seems unlikely this extra energy will remain in the 
form of pulses like those of figure 2. We know such pulses 
propagate unchanged in a flat space, but when they have 
traveled cosmological distances, so the curvature becomes 
noticeable, we should expect a gradual thcrmalization. 
We will model this thermalization by including a small, 
constant conductivity, a, in our equations. Note that a 
is not directly related to the normal conductivity of the 
plasma, which (in our cosmological units) would be very 
high, a is just a device to represent the slow thermaliza- 
tion, and will have a value of order unity. For technical 
reasons, which we explain below, we choose a = l/(2n). 

There are two reasons why we have to investigate the 
effect of a: 

1. The thcrmalization of the pulses raises the temper- 
ature of the plasma, and so tends to reduce the 
value of the ratio p. 

2. The slow decline in the height of the pulses will 
probably reduce the value of A in the tail, and this 
will also reduce the value of p. 

The gradual transfer of energy to the plasma implies 
that the product R(r])T(r]) will not remain constant, as 
in a simple expansion, but will slowly increase. In this 
respect a simulates inflation, but on a much longer time 
scale. 



XII. INCLUSION OF CONDUCTIVITY: FIELDS 

The wave number, k, the permittivity, e(n), and the 
conductivity, a, satisfy the dispersion relation 

k 2 = n 2 e(n) = n(n + Anicr) . (65) 

By analogy with Maxwell's equations in flat space, the 
equations (17), (18) and (19) for our dipole fields become: 

me sinh 2 x fi - 2/ 3 = 0, (66) 
inef 2 -U = 0, (67) 
fi - U + inh = 0. (68) 



A. Magnetic field 

From these we obtain an equation for / 3 alone: 

+ (" 2 '-^) /3 = a m 

Writing k 2 for n 2 e this equation takes a familiar form; 
f 3 can be derived from the empty-space formula sim- 
ply by writing k for n everywhere, except, of course, 
for the normalization function, which we denote now by 
C 3 (n,a): 

2u 

C 3 (n, a) is most easily determined in flat space; this is 
sufficient since we only have to consider small distances. 
The calculation is done in Appendix C; see (C18). 

B. Electric fields 

(66) and (70) give the radial component of the electric 
field: 

j _ -iC 3 (n,a) exp(^fcx)(l - u 2 ) 2 
4(n + Airi(j)u 3 

x (l - 2iku + u 2 ) . (71) 

(67) and (70) give the transverse component of the 
electric field: 

iC 3 (n, a) exp(zfcx) 

1 2 — 

4(n + 47t«(t)m 2 
x [(1 - w 2 ) 2 - 2iku{\ + u 2 ) - Ak 2 u 2 } . (72) 

C. Propagation of the magnetic field 

The Fourier transform of the magnetic field becomes 

— 3qV(n + 4iria)Ni (9) sin(nr) cxp(ifcx — inrj) 
12 (3n + 8ni<r)nu 

x (1 - 2iku + u 2 ) . (73) 

We can integrate around the pole at n = in the same 
way as before, except that we have to respect the branch 
points of k, at n — and n = —inia. We have also to 
take account of the pole at n = —8iri<r/3. A suitable 
contour is shown in figure 6. 

The integrals are straightforward, and the resulting 
propagation of F\ 2 is shown in figure 7. Notice that the 
pulses now show a tail that represents a reflected wave. 
This is unlikely to be significant in practice, because any 
irregularities in the plasma will tend to disrupt the co- 
herence of the wave as it converges on the origin. 

We note here for future reference that when r\ 3> 1 
we can derive a simple asymptotic form for F\ 2 and F20 , 
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FIG. 6: Contour for computing the magnetic field, F12, 
when conductivity is included. The points on the negative 
imaginary axis labeled A and B are at n = — 8nia/3 and 
n = — 4-Kia, respectively. The contour is traversed in a clock- 
wise sense. 



using the fact that in the integration around the contour 
the integrand is concentrated near n = 0. We just give 
the result: 



F\2 

F20 



iqV sin t 



rj + T 



exp 



rj + T 



F 12 /{Attu). 



(74) 
(75) 



XIII. INCLUSION OF CONDUCTIVITY: 
POTENTIALS 

Here we follow the prescription of Landau and Lifshitz 
[26] : use the same equations relating potentials and fields 
as in empty space, i.e. (39), (40) and (41), but modify 
the Lorenz condition by including the permittivity in the 
ho term: 



d 



znesinh 2 y/i + — Isinh 2 y/ii] — 2/i 2 
dx 



0. 



(76) 



A. Scalar potential 

From these equations, as before, we can derive an equa- 
tion for ho alone: 



d 2 h dh 
2- + 2cothx-^- + 



d X 2 



dx 



n 2 e 



sinh x 



h = 0. (77) 



This is the same equation as we obtained for empty space, 
except that, just as for the magnetic field, in place of n 2 
we must write k 2 = n 2 e. (49) applies as before, provided 
we write a = \Jk 2 — 1, and use a normalization function 
C±{n, a) in place of C2(n). 

This normalization function can be found from (40) 
and (71), in the limit x ^ 0. We get, in place of (50), 



din, a) 



iC 3 (n, a) 
4(n + Airier) ' 



(78) 




3 ^ 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 
X - 1 2 3 4 5 



For the Fourier synthesis we need, in general, to use a 
more complicated contour than the one shown in figure 6, 
because we have to take account of the branch points of 
both k and a. With our choice of a — l/2n, however, the 
branch points of a coalesce into a single point at n = —i, 
and the contour of figure 6 is adequate. 



B. Vector potential 

(52) holds, as in the case of zero conductivity. The 
first term on the right side of (52) comes from (72): 



-if 2 _ C 3 (n,<r)exp(ikx) 



x [(1 - u 2 ) 2 - 2iku(l + u 2 ) - Ak 2 u 2 } . (79) 

The second term on the right side of (52) comes from 
(49): 



FIG. 7: Propagation of the magnetic field, F12, in the presence 
of conductivity, a — l/2n. 



iho iCi(n, a) exp(iax)(l — u 2 ) 
n nu 2 

x (l - 2iau + u 2 ) . (80) 
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The two normalization functions arc related by (78), 
and we also know k 2 = n(n + Anicr). Combining, we get: 



h 2 



C 3 (n, a) 



(81) 



4n(n + 4iria)u 2 
K a = e lk * [(1 - u 2 ) 2 - 2iku(l + u 2 ) - Ak 2 u 2 } , (82) 
A a = e lax (l-u 2 )(-l + 2iau-u 2 ) . (83) 



D. Asymptotic form of the vector potential 

We can show easily that in the limit i] — > oo, A 2 sat- 
isfies a diffusion equation. The main terms come from 
the integral down the cut, and in this limit the integrand 
is concentrated close to n = 0. This makes it possible, 
as for the fields, to obtain a simple expression for the 
asymptotic form: 



C. Propagation of the vector potential 

Combining (81) with (37), (C18) and (10), and setting 
— 87ricr/3 = n p , we get 

^ = -,m Wsi n( T ) r ^ 

2(n — n p )n z u z 

We note that going from (53) to (84) takes three simple 
steps: 

1. Rename N and A] call them JC a and A a , respec- 
tively. 

2. In the final parenthesis, substitute k for n i n both 
/C CT and Aa ; t his inclu des redefining a to be \Jk 2 — 1 
rather than \/n 2 — 1. 

3. In the prefactor, change n 3 in the denominator to 
n 2 (n — n p ). 

(84) becomes (53) in the limit a — > 0, as it should. 

We transform to rj, x, 9, </> coordinates as before, by 
dividing (84) by 2tt and integrating clockwise around the 
contour of figure 6. The result is shown in figure 8. 



3qVr sin6 
A 2 « F(a,x,v), 




i i i i i i i i i i i i i i i i i i i i i i i i i i i i 

X - 1 2 3 4 5 6 

FIG. 8: Propagation of the vector potential, A2, in the pres- 
ence of conductivity, a = l/2n. 



16tt 

(l-u 2 ) 2 



(85) 



erf I X\ l — ' - 



(tu^ \ ' " y 77 / a Y l 2 r\ z l 2 u 
x [(1 + m 2 )?? + 47TcrMx] exp ( ~™ X ^ } . (86) 



/ J 

When x is also large, of order ^/rj, the term involving 
the error function is negligible compared to the others. 

In figure 9 we plot ^JrjA 2 against £ = x/^/Vj with A 2 
computed using the asymptotic formula. We see y^rjA 2 
tending to a constant form; with our choice of abscissa, 
this is a Gaussian for large \ and 77. For 77 > 4, this 
graph is indistinguishable from the one using integration 
around the contour. 




FIG. 9: yfrjAi plotted against £ = x/\^j to demonstrate the 
asymptotic form. In this graph we use the asymptotic form 
(85) for the tail. 



XIV. THE ENERGY-TRANSFER EQUATION 

We will now investigate the effect of a non-zero con- 
ductivity on the temperature of the plasma. We start 
from Weinberg [6], equation (4.7.9): 

Tr„ =-LJL {s/gT^) + r^T^. (87) 
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This expression is zero for any value of v, but the most 
important one for us is v = 0. We treat the plasma as 
a perfect fluid at rest in FRW coordinates rj, \, 9, 0, so 
that the 4-velocity (U°, W) has U l = 0, i = 1, 2, 3. Since 
g^UW = -1 and .goo - R 2 (v), U° = 1/R. 

Multiply (87) by ^/gd A xU Ul sum over v, and set the 
resulting expression equal to zero; we get 



^gd 4 xU T^ = 0. 



(88) 



The left-hand side of this equation is a coordinate scalar 
with ^/gd 4 x an element of proper volume. Expanding 
(88): 



y/gd 4 xU 



1 



d 



+ r%T ij + t° 00 t 00 } = o 



(89) 



From this we can derive an integrodifferential equa- 
tion for T 00 . Details are given in appendix D. The 
resulting equation is most simply written in terms of 
Am(v) = R 4 {v)T 4 (r]) (eV) 4 • m 4 . A similar quantity was 
introduced in section X, but there it could be treated sim- 
ply as a constant. Here we calculate its variation with 
time: 



, a(rj - rj') 



2+p( V ) 



x exp [— 27177(77 — v')] = 0, 
P 4 = 3.97x lO^cV-m. 



(90) 
(91) 



Note that because of the appearance of the function 
p(rj), this equation by itself is not sufficient to calcu- 
late A M . We obtain an independent equation connecting 
these two functions in the next section. 



XV. MASS GENERATION: SECOND 
CALCULATION 

We will now obtain an equation for the ratio p(r]) anal- 
ogous to (64), but based on the asymptotic form (85) for 
the tail of the vector potential. In our previous calcula- 
tion, where the tail of A 2 reached a constant asymptotic 
value, we could be sure that p would rise monotonically; 
the only question was whether there would be factors of 
Am that would make the increase far too fast or far too 
slow. (This is just what happens, for example, if we use 
quadrupole rather than dipole potentials.) In this second 
calculation, we can anticipate that there will not be any 
unwanted factors of Am, but only a detailed calculation 
will show whether p continues to rise with 77. 

The analog of (59) is 



diV s = 47r(0.24)4 1 /4 (7 ?s ) sinh 2 X d X . 



(92) 



To get the contribution to A\ at rj p from the particles 
in the shell, we multiply (92) by A\\a, x, (r] p — r] s )], with 
A 2(^,X,V) given by (85). 



Proceeding as before we get the analog of (62): 
2*(0.24)£ p / < 4 fa) \ 

7p ) j Arh [2+p(rk) r 



3i? 2 (r? p 



x / d X .F [<t,X,(Vp-Vb)], 
Jo 

with Fo given by (86). 

The collision time, t, will be given by 

(l+p) 



(93) 



87rln A(0.24)g 2 < 4 (77 s )(2 +p)3/ 2 ' 



(94) 



We can now get an expression for p by using (94) in 
(93) and dividing by T 2 (?7p): 

o 

p(v P ) 



2567r 2 lnA4f (r/ p ) 



M VIPJ 



x / d X F 2 [a,x,(r? P -r, s )], (95) 
Jo 

with F given by (86). 

A. Calculation of p and Am 

As before, we set In A = 10 for dcfinitcness. Simultane- 
ous numerical solution of (90) and (95) is straightforward, 
but we need to choose a starting value for Am- We have 
used two extreme values: 

1. A M starts at 2.5 x 10~ 26 eV 4 • m 4 . This is a "nat- 
ural" value, in the sense that the Universe begins 
with the typical wavelength of the radiation ap- 
proximately equal to the radius, R. p becomes 
equal to unity around i] cw = 1650. 

2. A M starts at 2.5 x 10 74 eV 4 • m 4 . This is approx- 
imately equal to the observed value today, p be- 
comes equal to unity around rj cw — 480. 

These values of rj cw are significantly larger than the 
previous value of 11 when a = 0, but they are still within 
a few orders of magnitude of unity. The most important 
feature of this calculation is that it shows p(rj) does con- 
tinue to rise, even when a ^ 0, so a CW transition will 
still eventually become possible. 

XVI. CONCLUSION 

We have considered the problem of the generation of 
a particle mass scale in a theory in which this scale is 
not determined by the gravitational constant, through 
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Wpianck- Conformal gravity [2] may turn out to be a the- 
ory of this type. We find that in such a theory distant 
matter is significant in an unexpected way and the mass 
scale is to a large extent classically determined. Since 
the mass scale develops through the agency of the classi- 
cal vector potential, it is bound to appear in any model 
universe containing charged particles, providing the un- 
derlying space has sufficient curvature and allows enough 
time. 

Many questions remain, among them the following: 

1. Will a CW transition necessarily take place as the 
ratio (A 2 ) nt /(A 2 ) th approaches unity? 

2. Is the tail of A important only in the early Uni- 
verse, or are there other places, regions of high 
gravitational fields, where its effects can be ob- 
served even now? 

3. How does the real complicated early plasma deter- 
mine such things as the collision time? 

4. Finally, perhaps most important, can we really ex- 
pect A to propagate over cosmological distances? 
We have taken findings from ordinary plasma the- 
ory and used them in the very different circum- 
stances of the early Universe. 

These considerations are beyond the scope of this pa- 
per, which is solely concerned with the interplay, in a 
simple model, of quantum mechanics (density matrix, 
the Lagrangian for a scalar field, the Coleman- Weinberg 
transition) and the classical equations of propagation of 
the ordinary vector potential in a FRW space of negative 
curvature. 
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APPENDIX A: HYPERGEOMETRIC 
FUNCTIONS 

Hypergeometric functions are usually difficult to han- 
dle because they depend on three parameters. For all the 
hypergeometric functions we need in this paper, however, 
the parameters a, b, and c can be expressed in terms of 
two numbers, an integer, d, and a quantity a which is 
complex in general: 

a = d + ice; b = d — ice; c = d+ 1/2. (Al) 



Hypergeometric functions that satisfy this condition 
can be expressed in terms of elementary functions. 
For proof we use formulae from Abramowitz and Ste- 
gun ([16], hereafter AS). Define fi = (1 — cosh(x))/2. 
Then, for d = 0, AS (15.1.17) gives, in our nota- 
tion, F(ia, — ice; 1/2; \i) = F[ice, — ice; 1/2; sin 2 (ix/2)] = 
cos(ax). Other functions can be obtained by successive 
differentiation or integration using AS (15.2.1). If a hap- 
pens to be purely real or purely imaginary the hyperge- 
ometric functions are real. 

For dipole propagation the only functions we need are 
those with d = — 1 and d = 2. They are listed below, 
together with the functions with intermediate values of 
d. We can replace the argument fi by \: 

F(-l + ia,-l- ice; -1/2; X ) = 
cosh x cos(cex) + ct sinh x sin^x) > ( A2) 

F(ia,-ia;l/2; X ) = cos(ax), (A3) 

F(l + ia,l-ia;3/2; X ) = ^4^, (A4) 

a sinh x 

F(2 + ia,2-ia;5/2; X ) = 
— 3[asinhxcos(ax) — coshx sin(ax)] 



a(l + a 2 ) sinh x 



(A5) 



APPENDIX B: PARAMETERS IN THE MODEL 
OF MANNHEIM 

Mannheim finds that the field equations of conformal 
gravity take on a particularly simple form in the context 
of a FRW space, as in cosmological investigations. He 
obtains ([2], equation (230)) the following expression for 
the expansion factor, i?, as a function of the ordinary 
time, t, assuming the matter in the Universe is in the 
form of radiation and the space has negative curvature 
(fc = -l): 



R 2 {t) 



(/3 M - 1) /3 M sinh 2 (a;^ 2 rf) 



2otM 



(Bl) 



«m and /3m are real, positive parameters. /3m is di- 
mensionless and greater than unity; «m has dimension 
length -2 . (We write cum rather than Mannheim's a to 
prevent any possible confusion if a is used elsewhere, for 
example to denote the fine structure constant; we use 
the same convention for other parameters in Mannheim's 
model.) We will use (Bl) for the whole range of time from 
t = to the present, recognizing that this will be inaccu- 
rate during the mass-dominated phase of the expansion. 
We estimate this time interval below. 

Time t = corresponds to the surface of the bub- 
ble. Here R has a non-zero minimum value, i? m ; n = 
[(/3m — 1)/(2o;m)] 1 ^ 2 , and a maximum temperature, T max . 
Mannheim defines a third parameter, Am, by 



A M = R\t)T A (t). 



(B2) 



We will not try to explain the large size of Am , but will 
make essential use of the fact that (in this simple model, 
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at least) it has a constant value during the expansion. If 
we can determine the three parameters om, /3m and 
we can derive the values of all the rest of Mannheim's 
parameters for this model. 

Equation (Bl) leads to an expression for the acceler- 
ation parameter, qo. By comparing this with the value 
obtained from supernova observations, Mannheim finds 
that a^ct k 0.7 at the present time, so that a\j^ is of 
the order of the inverse Hubble distance. Am can be esti- 
mated from the current CMB spectrum and the Hubble 
distance; we will use a value of A^ 4 = 1.68 x 10 29 . 

The parameter (3m is harder to estimate. The temper- 
ature evolution of the Universe is given by the following 
equation ([2], equation (232)): 



rp2 



= 1 + 



2f3 M sinh 2 {a^ct) 
(Pm - 1) 



(B3) 



Since T max S> T , where T is the current temperature, 
(3m must be extremely close to 1. It is therefore con- 
venient to define a new parameter, i5m = (3m ~ 1> with 
< <5m <C 1. Equation (Bl) then shows that 6m has an 
appreciable effect on R(t) only at the very earliest times. 
We will estimate 5m below. 

Mannheim suggests that the CW phase transition oc- 
curs at an intermediate temperature, TV, with T v « 
T max T . He regards the mechanism of this phase transi- 
tion as a matter for particle physics, and does not con- 
sider it in detail. In our model the transition occurs at a 
considerably higher temperature, and appears not to be 
connected to Ty. 

Define x = a 1,>2 t, with xq = a 1 / 2 ^ = 0.7 the cur- 
rent value. The current CMB temperature is To = 
2.3 x 10~ 4 eV. The CW transition will occur at a tem- 
perature of about T cw = 10 9 eV, so, since 8m <C x cw <C 1: 



Vo - lew 



x cw = sinh(0.7)(2.3 x 10~ 4 )/10 9 
= 1.7 x 10~ 13 , 
x ° dx 



-L 



„ sinMx) 
h, I tanh (*o/2) . = 2{) _ 



Also, since x cw <C 1, 

Vc 



tanh(x cw /2) 



da: 



(B4) 



(B5) 



o [S M /2 + x 2 } 1/2 



= In 



x + + x- 

V^m/2 

In (2 3 / 2 a;cw( 5 M 1/2 ) , 
2- 3 / 2 < 2 cxp(, 7cw ), 

0-1/2x1/2 
Z °M ' 



(B6) 

(B7) 
(B8) 



^cw^cw/ *£min 

0.5 x 10 9 x cxp(77 cw ) = 3 x 10 13 eV, (B9) 



where we have used r/ cv/ — 11, as derived from (64). 
Equation (232) of [2] gives 

S M - 2(0.76) 2 ( ™*™~ 4 ) = 6.8 x 10" 35 , (BIO) 



(Bll) 



v 3 x 10 13 J 
T v = T max (<5 M /2) 1/4 = 7 x 10 4 cV 



Our equation (Bl), Mannheim's expression for R 2 (t), 
will become inaccurate during the mass-dominated phase 
of the expansion, between T = T cw and Ty. This cor- 
responds to an interval of conformal time of about 10, 
compared to 40 for the whole interval from t — to the 
present. 



APPENDIX C: NORMALIZATION WITH 
CONDUCTIVITY 

In this appendix we will work in the usual spherical po- 
lar coordinates, and make connection with Ricmannian 
coordinates when necessary. 

Suppose we have a dipole at the origin, oscillating with 
time dependence exp(— iut) in the z direction. The sur- 
rounding medium is of uniform conductivity, cf, so that 
the current density, j, is given by j = <rE. We will ana- 
lyze this system by imagining a small sphere of radius n 
cut out of the medium surrounding the dipole. Induced 
currents flowing in the medium will cause surface charges 
to appear on the sphere, and the total dipole moment will 
be the sum of the original dipole moment and that due 
to the induced charges. We assume the permittivity and 
magnetic susceptibility are essentially unity, so D = E 
and B = H. In such a system V • j = follows from 
Maxwell's equations, so there are no volume charges in 
the medium. 

The dipole moment at the center of the small sphere is 
denoted by D true — D(uj) exp(-iu)t), where D(w) is the 
true dipole strength at angular frequency u>. The induced 
dipole moment due to the surface charges is D- m d, so the 
total dipole moment is D tot = D true + And- 

Just outside the sphere the electrostatic potential is 
given by 



$ = D tot Pi (cos 6) /rl, 



(CI) 



where Pi (cos 9) — cosO. The radial component of E is 
given by 

E r = 2 D tot P 1 (cos 8)/r\. (C2) 
The surface charge density, s, obeys the relation 
ds 



dt 



= -Jr 



(C3) 



where j r is evaluated just outside the sphere. This gives 
— iaE r 



s = 



to 



= s Pi (cos 9), 



(C4) 
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where 



-2iaD tot 
S ° = IE 



(C5) 



The induced dipole moment, And, is then given by an 
integral over the surface of the sphere: 



P7T 

D ind = 2ir / e*0 sin 0r 2 [s o A (cos 60][r? Pi (cos 0)] 
Jo 



-Airia\ ( 2 



giving 



Dtot — Artie + And 

3lo + 8iria 



(C6) 



(C7) 



We set i/^ = C 5 (w, a)Ni(6)h(r) exp(— iwi), where 
C 5 (cu,a) is the normalizing function we are looking for, 
Ni = dPi/d6, and h{r) satisfies 

d 2 h 2dh 2h 2 , „ 

— + -- 2 - + K 2 / i -0, C8 

dr z r dr r z 

with k 2 = w(w + Airia). We choose the solution that 
represents outgoing waves, so 



h{r) = h ( i\Kr 



(fir) 2 nr 



(C9) 



Here /i^(/tr) is the spherical Bessel function defined in 
AS, ch. 10. 

The Maxwell equation VxH = + AniajE then 
gives 

E = -2iC 5 (oj,W)Pi (cos e)h[ 1] (kt) expj-iut) 
r(tu + Ania) 

For small r this becomes: 

— 265(0;, cr)Pi(cos0) exp(— iwi) 



E r = 

K 2 r 3 (uj + Airia) 

But also, for small r, we have (C2), so 

-3wk 2 (w + Airia)D{uj) 



H<p = 
E r = 



(Cll) 



(C12) 



(3w + 87ria) 
-3wk 2 (w + 47rio 7 ) £>(^) 
(3w + 87riCT) 

xAi(6»)/i* 1) (Kr)exp(-iwt), (C13) 



r(3w + Siria) 



(C14) 



In this appendix we have used ordinary polar coordi- 
nates in flat space. We need now to transform to coordi- 
nates of a FRW flat space, with metric 



ds 2 = i? 2 (7 7 )(-dr / 2 +d X 2 +X 2 d0 2 



+ X 2 sin 2 6»d0 2 ) 



(C15) 



To get from H4 of (C13) to F 12 of (13), in the flat 
FRW metric (CI 5), we first follow Weinberg [6], ch. 4, 
sec. 8, and multiply by \. We then convert from t to r) 
and r to \ by multiplying by R 2 (rj). We also convert w 
to n, k to k and a to a = R(r})a: 



'12 



-3nfc 2 (n + Awia)D(n) 



R(r])(3n + 8iTia) 

xN 1 {e)xh[ 1) (k X )eM~inv)- (C16) 

We can obtain the corresponding function in curved 
space by multiplying (70) by Ni(6) exp(— inrj): 
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C7 3 (n,a)JVi(g)exp[i(fcx-n77)] 
2w 

x (l - 2ifcu + u 2 ) . 



(C17) 



We can now find Cs(n, a) by matching (C16) with 
(C17) for small X - 



C 3 (n,a) = 



3m(n + Ama)D{n) 
R(r])(3n + 8iria) 



(C18) 



With D(n) given by (10), we can check that in the 
limit cr -> 0, (C18) tends to (32), as it should. 



APPENDIX D: DERIVATION OF THE 
ENERGY-TRANSFER EQUATION 

In this appendix we present details of the derivation of 
the energy-transfer equation (90), starting from equation 
(89): 



y/g d 4 x Uq 



+ r°,T" + r° t 00 ] 



0. 



Write this as 



z\ + z 2 + z 3 + z 4 = 0, 



(Dl) 



(D2) 



where the are derived from the four terms in the 
bracket in (Dl). The strategy will be to get an expres- 
sion for Z2 from the propagation of E and H. We combine 
this with z\, Z3 and Z4, and arrive at an intcgrodifferen- 
tial equation for T 00 . 



1. Calculation of 21 

We consider first the effect of a single source particle 
at the origin of coordinates, and let the duration of our 
4- volume be Aa: . Our coordinate volume dxd9d(j) is a 
small part of a spherical shell of radius \ surrounding the 
source particle. Then z\ in (D2) is given by 

Zl = Ax°U -^(y/gT 00 )dxd9d<j> 

= P 4 (ART 00 + PT 00 ) sinh 2 X sin eArjdxded^. (D3) 
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2. Calculation of 2 2 

For z-i in (D2) the only relevant derivative is d/d X , and 
we get 

z 2 = Ax d9d4>U A(y^T w ) 

= Aryi?— (Vff? 110 ) d X d9d<f>. (D4) 

Here A(y / <?T 10 ) is the net outward flux from the shell, 
for a pulse that protrudes on both sides. The pulse is, 
of course, attenuated on the outward side, so the total 
outward flux is negative. 



3. Calculation of z 3 and za 

To compute z 3 and Z4 in (D2) we start from Weinberg 
(5.4.2): 

T"" = pg^ + (p + p)U»U v . (D5) 

The only Christoffel symbols we need are Tq = R/R 
and T° = Rgij/R 3 , where a dot denotes d/dr], and i, j 
take values 1, 2, 3. Then 



T oo zE 1 p + p = JL 

R 2 i? 2 i? 2 ' 

T ij = pg ij , 

r o ™a Rp Rpgij9 ij 
^ A r 3 i? 3 



2pR 



R 3 



if p = p/3. 



(D6) 
(D7) 

(D8) 



We now have 



2 pR. 



Z3 + Z4 = U{)—j^-Axdx&Od(j)^/g 
= 2T oa RAT 1 ^d X d9d(f). 



4. Energy-transfer equation in terms of T p 



Using (D3), (D4) and (D9) in (D2) we get 



(D9) 



4i?T 00 + RT 00 + 



R 9 '-10^2 



sinh x ®X 



(T 10 sinh 2 X ) 



+ 2i?T 00 = 



^+T°° + -^^(T 10 sinh 2 x) = O.(DIO) 

We must now develop this equation in two ways: 

1. Express the first two terms as functions of the tem- 
perature, T{rj). 

2. Sum over all the source charges that contribute to 
T 10 . Since we are concerned here with the main 
pulse, not the tail, this implies an integration along 
the light cone. 

5. Expression of T 00 in terms of T(r)) 

We will measure temperature in electron-volts, R in 
meters, and will assume the energy density is given in 



terms of the temperature by the familiar formulae for 
electromagnetic radiation (we will clearly miss a numer- 
ical factor here when dealing with the real plasma, but 
the overall relationship will be correct). Denote by </mks 
the elementary charge in MKS units, and let the energy 
density be p. Then 



9 
Pi 



P1T 4 J • nT 



7T 

15 



tt 2 (1.6x 10- 19 )' 



15 (1.05 x 10 - 34 ) cl (3x 10 8 )' 



= 13.8 (eV) -J-nr 



rpOO 



J • m 



4P 1 T 3 (r ? )T(r ? ) 2P 1 T i ( V )R(r ] ) 



R*( V ) 



R 3 ( V ) 



(Dll) 
(D12) 
(D13) 

J-m- 5 .(D14) 



6. Integrating along the light cone 

Using the asymptotic forms (74) and (75) for the fields, 
we obtain 



-.10 



~9 9 9 ^12-^20 
(3gUsin6») 2 1 
i? 6 (r/) sinh 2 X 47I "(»7 + T ) 

A(rio s inh 2 X ) = 



exp 



,(D15) 



-('SqV sm9) 2 ax 
R 6 (V)(V + T) 2 



■ exp 



-2nar 

T] + T 



.(D16) 



This form refers to a single source particle. To get the 
total T 10 term in (D10) we must integrate along the light 
cone. We will need the number density, which we write as 
P2T 3 , with the constant P 2 defined as for electromagnetic 
radiation: 



ft = Ac(3)#f 



= 0.24 



1.6 x 10" 19 



(1.05 x IO- 34 ) (3 x 10 8 ) 

^20 



3.15 x 10 2U (eV) -nT 



(D17) 



We can now use (D16) to write an expression for the 
total T 10 term in (D10). We neglect r in comparison with 
77, and write 77 — t]' in place of r\. We also express the 
charge q, which is in Gaussian units, in terms of (?mks : 



\/47reo 



<?mks ■ 



(D18) 



with e = (1/36tt) x 1Q- 9 C 2 • J" 1 
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d 



, (T 10 sinh 2 X ) 
snhr x °X 

x (47ri? 3 (?7')sinh 2 x) 



i? 6 (r/) 



-, pr) noo 



sinh 2 x 



(3 9M Ks) 2 (^ 2 )(si^ 9)ax 



47re (?7 - f?') 2 



exp 



-27nr X 2 
77-77' 



Ar?. 



(D19) 



As in section X we set (V 2 ) w 1/(2 +p(r/)) and (sin 2 



2/3; we also define a new constant P 3 : 



F 3 = P 2 x 36tt x (2/3) x (1.6 x 1(T 19 ) 2 x (9 x 10 9 ) 



5.47 x 1(T 6 (oV) -3 • J • m~ 2 , 



d 



sinh 2 x d X 



T 10 sinh 2 X ) 



■Pa 



R 6 (v) Jo 



dr] 



f T 3 ( V ')R 3 ( V ') 
2+p{r]) 



a(r] — t] ') exp [— 2ncr(rj — 77 ')] A?7 . 



(D20) 
(D21) 



r 



7. Final form of the energy-transfer equation 

Using (D13), (D14) and (D21) in (D10): 

4P 1 T 3 (r/)T(?7)i? 4 (77) + 4P 1 T 4 (r ? )i? 3 (r/)i?(r/) 
r\ ,T 3 (r?')iW) 



able Am(v) = T 4 (ri)R 4 (r]) eV 4 • m 4 , and a new constant, 
P 4 : 



A M {v) 



P 3 



/ dry': 
Jo 



2+p(rj) 
x exp [-27rcr(?7 — n')\ = 0. 



•(7(77-77') 



Jo 



(7?) 2 + ^) 



(D22) 



x exp [-27ra(77 - rj ')] = , (D23) 
P 4 = P3/-P1 = 3.97 x 10" 7 eV • m . (D24) 



This equation can be written in terms of the single vari- 
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